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The presence of long-range quantum spin correlations un- 
derlies a variety of physical phenomena in condensed mat- 
ter systems, potentially including high-temperature supercon- 
ductivity 1 1 , 2|. However, many properties of exotic strongly 
correlated spin systems (e.g., spin liquids) have proved diffi- 
cult to study, in part because calculations involving A/^-body 
entanglement become intractable for as few as A/^ ~ 30 par- 
ticles [3J. Feynman divined that a quantum simulator — a 
special-purpose "analog" processor built using quantum par- 
ticles (qubits) — would be inherently adept at such problems 
14, 5 1. In the context of quantum magnetism, a number of ex- 
periments have demonstrated the feasibility of this approach 
ISHB I . However, simulations of quantum magnetism allow- 
ing controlled, tunable interactions between spins localized 
on 2D and 3D lattices of more than a few lO's of qubits have 
yet to be demonstrated, owing in part to the technical chal- 
lenge of realizing large-scale qubit arrays. Here we demon- 
strate a variable-range Ising-type spin-spin interaction Jij on 
a naturally occurring 2D triangular crystal lattice of hundreds 
of spin- 1/2 particles (^Be+ ions stored in a Penning trap), a 
computationally relevant scale more than an order of magni- 
tude larger than existing experiments. We show that a spin- 
dependent optical dipole force can produce an antiferromag- 
netic interaction Jij oc jrj^, where a is tunable over < a < 3; 
dij is the distance between spin pairs. These power-laws cor- 
respond physically to infinite-range {a = 0), Coulomb-like 
(a = 1), monopole-dipole (a = 2) and dipole-dipole (a = 3) 
couplings. Experimentally, we demonstrate excellent agree- 
ment with theory for 0.05 <a<lA. This demonstration cou- 
pled with the high spin-count, excellent quantum control and 
low technical complexity of the Penning trap brings within 
reach simulation of interesting and otherwise computationally 
intractable problems in quantum magnetism. 
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Figure 1. The Penning trap confines hundreds of spin- 1/2 particles 
(qubits) on a two-dimensional (2D) triangular lattice. Each qubit is 
the valence electron spin ofa9Be+ ion. (lower) A Penning trap con- 
fines ions by use of a combination of static electric and magnetic 
fields. The trap parameters are configured so that laser-cooled ions 
form a triangular 2D crystal. A general spin-spin interaction Hj is 
generated by a spin-dependent excitation of the transverse (along z) 
motional modes of the ion crystal. This coupling is implemented 
with an optical dipole force (ODF) due to a pair of off-resonance 
laser beams (left side) with angular separation Qr and difference fre- 
quency jUR. Microwaves at 124 GHz are directed to the ions with a 
horn and permit global spin rotations Hb. (upper) A representative 
top-view resonance-fluorescence image showing the center region of 
an ion crystal captured in the ions' rest-frame; in the lab-frame the 
ions rotate a at (Or= 27i x 43.8 kHz |21 1. Fluorescence is an indica- 
tion of the qubit spin-state (|t) bright, ||) dark); here, the ions are in 
the state |t). The lattice constant is do ~ 20 /jm. 



A challenge in condensed matter physics is the fact that 
many quantum-magnetic interactions cannot currently be 
modeled in a meaningful way. A canonical example is the 
spin-liquid postulated by Anderson 1 1 1. He suggested this ex- 
otic state would arise in a collection of spin- 1/2 particles resid- 
ing on a triangular lattice and coupled by a nearest-neighbor 
antiferromagnetic Heisenberg interaction. The spin-liquid's 



ground state is massively degenerate, owing to spin frustra- 
tion, and is expected to exhibit unusual behaviors including 
phase transitions at zero temperature driven by quantum fluc- 
tuations ifTSl . However, despite recent advances (161 1121 a 
detailed understanding of large-scale frustration in solids re- 
mains elusive [.2l [T8H^ . 
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Atomic physicists have recently entered the fray, providing 
a bottom-up approach by engineering the relevant spin inter- 
actions in quantum simulators |[5]|22l[23l. The necessary ex- 
perimental capabilities were first demonstrated in the context 
of atomic clocks: laser cooling, deterministic spin localiza- 
tion, precise spin-state quantum control, high-fidelity readout, 
and engineered spin-spin coupling (e.g., |24|). In the domain 
of quantum magnetism, this tool set permits control of param- 
eters commonly viewed as immutable in natural solid systems, 
e.g., lattice spacing and geometry, and spin-spin interaction 
strength and range. 

Initial simulations of quantum Ising and Heisenberg inter- 
actions with localized spins were done with neutral atoms in 
optical lattices ElEl, atomic ions in Paul traps ll9l fTQlfT3l[T4l . 
and photons \Vl]. This work involved a domain readily cal- 
culable on a classical computer: interactions between ~ 10 
qubits localized in ID chains. The move to quantum magnetic 
interactions on 2D lattices and larger, computationally rele- 
vant particle numbers is the crucial next step but can require 
significant technological development f25\ . 

In our Penning trap apparatus, laser-cooled ^Be+ ions natu- 
rally form a stable 2D Coulomb crystal on a triangular lattice 
with ~ 300 spins (Fig.[T]). Each ion is spin- 1/2 system (qubit) 
over which we exert high fidelity quantum control |26|. In 
this paper we demonstrate the use of a spin-dependent optical 
dipole force (ODF) to engineer a continuously-tunable Ising- 
type spin- spin coupling Jij oc jrj^. This capability, in tandem 
with a modified measurement routine (e.g., by more sophisti- 
cated processing of images like Fig. [T]), is a key advance to- 
ward useful simulations of quantum magnetism. 

A Penning trap confines ions in a static quadrupolar electric 
potential (see Methods) and a strong, homogeneous magnetic 
field Bqz (here, = 4.46 Tesla). Axial trapping (along z) is 
due to the electric field. Ion rotation at frequency co^ (about 
z) produces a radial restoring potential due to the velocity- 
dependent Lorentz force (qv x B). Tuning the ratio of the ax- 
ial to radial confinement allows controlled formation of a pla- 
nar geometry and, after Doppler laser cooling, the formation 
of a 2D Coulomb crystal on a triangular lattice (see 1211 and 
Methods). We routinely generate crystals with 100 < A/^ < 300 
ions, where the valence-electron spin state of each ion serves 
as a qubit |26|. Following techniques developed in linear (ID) 
Paul traps |27|, spins confined in the same trapping potential 
are coupled through their shared motional degrees of freedom. 

Using well controlled external fields, we engineer spin in- 
teractions of the form 



(1) 



where is the z-Pauli matrix for ion /. We label the qubit spin 
states It) \ms = +1/2) and ||) = -1/2), where is 

the spin's projection along the Bqz quantizing field, so that 
It-) = + It-) and \U) = - ||,-). Hb is an interaction (e.g., 
B/jCx) generated by externally applied microwaves at 124 GHz 
that couples equally to all spins and permits global rotations 
(Fig.[T|). The interaction Hi describes a general coupling Jij 
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Figure 2. Spin-spin interactions are mediated by the ion crystal's 
transverse motional degrees of freedom, (a) For a 2D crystal with 
N = 2\1 ions and cOr = 271 x 45.6 kHz, we calculate the eigenfunc- 
tions byn and eigenfrequencies co^ for the N transverse motional 
modes (see Supplementary Information). Plotted here are oOm and 
bm for the 14 highest frequency modes. Relative mode amplitude is 
indicated by color. The center-of-mass motion (COM) is the highest 
in frequency (coi ~ 27i x 795 kHz) and has no spatial variation in b^- 
The lowest-frequency mode (CO217 ~ 27i x 200 kHz) exhibits spatial 
variation in bm at the lattice spacing length scale ^ 20;wm. (b) By 
use of Eq.jS] Jij is calculated explicitly for A/^ = 217 spins and plotted 
as a function of spin-spin separation dij. For — coi < 271 x 1 kHz, 
Hqdf principally excites COM motion in which all ions equally par- 
ticipate: the spin- spin interaction is spatially uniform. As the detun- 
ing is increased, modes of higher spatial frequencies participate in the 
interaction and Jij develops a finite interaction length. We find the 
scaling of Jij with spin separation dij follows the power law d^^. 
For fdR — QOi 2nx 500kHz, all transverse modes participate and 
the spin- spin coupling power law approaches a = 3. The solid lines 
are power-law fits to the theory points. For comparison with other 
experiments, the nearest neighbor coupling (do = 20 /jm) is denoted 
by the dashed line, (c)-(e) The power-law nature of Jij is qualita- 
tively illustrated forN= 19 (larger N is illegible). Nodes (spins) are 
joined by lines colored in proportion to their coupling strength for 
various a. (f) For context, the graph for a ID nearest-neighbor Ising 
interaction, a well-known model in quantum field theory. 
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between spins / and j, distance dfj apart ||28l [29l ; Jij > 
(Jij < ) is an antiferromagnetic (ferromagnetic) coupling. 

We implement Hj using a spatially uniform spin-dependent 
optical dipole force (ODF) generated by a pair of off-resonant 
laser beams with difference frequency jur (see Fig.[T]and Sup- 
plementary Information). The ODF couples each ion's spin 
to one or more of the N transverse (along z) motional modes 
of the Coulomb crystal by forcing coherent displacements of 
the ions that in turn modify their Coulomb potential energy 
through the interaction 

N 

HoDF = -Y,F,{t)zial (2) 

i 

Here Fz{t) — Fq cos jURt is the ODF, Zi = 
I.^=ihm^/^ {ame-'''-' ^aie'''^') is the axial posi- 
tion operator for ion /, and bi^m are elements of the N 
transverse eigenfunctions bm at frequencies co^ normalized as 
1.^=1 Hm? = Iti \hm? = 1 12811291. The modes include 
the center-of-mass mode (coi) as well as an array of higher 
spatial-frequency modes that may be derived from atomistic 
calculations (see Fig. [2]) and confirmed by experimental 
measurement |30|. 

For small coherent displacements, where residual spin- 
motion entanglement can be neglected (see | 29 | and Meth- 
ods), HoDF is equivalent to Hi in Eq.[l] spins / and j are cou- 
pled in proportion to their spin states and d^- and mutual 
participation in each motional mode m. The coupling coeffi- 
cient is given by 1291 

These pairwise interaction coefficients Jij can be calculated 
explicitly by use of Eq. [3] and classical calculations of ion 
motional modes. We find that the range of interaction can 
be modified by detuning away from the COM mode as dis- 
cussed in Fig. [2] In the Hmit jur-o^i > 2;: x 500 kHz, all 
modes participate equally in the interaction and Jij oc d^j , as 
discussed in 1281 . At intermediate detuning, we find a power- 
law scaling of the interaction range Jij oc df^, where a can 
be tuned < a <3. That is, by adjusting the single experi- 
mental parameter jur we can mimic a continuum of physical 
couplings including important special cases: a = is infinite- 
range, a = 1 is monopole-monopole (Coulomb-like), a = 2 is 
monopole-dipole and a = 3 is dipole-dipole. Note that a = 
results in the so-called J^ interaction that gives rise to spin- 
squeezing and is used in quantum logic gates (see Supple- 
mentary Information) {TT\ . In addition, tuning jur also permits 
access to both antiferromagnetic (AFM, jur > o^i) and ferro- 
magnetic (FM, CO2 <C /i/? < COi ) couplings 1 13 |. 

Experimentally, we demonstrate a tunable-range Ising in- 
teraction by observing a global spin precession under appli- 
cation of Hi. We compare experimental data with the mean 
field (MF) prediction that the influence of Hj on spin j can be 
modeled as an excess magnetic field in the z-direction Bj due 
to the remaining N —1 spins, Bj = ^ Y^u^jJij (6^) (see Sup- 
plementary Information). For a general qubit superposition 
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Figure 3. Benchmarking the 2D Ising interaction, (a) Spin- 
precession benchmarking sequence for Hj. The spins are prepared 
at the outset in |t) (a ferromagnetic state). The spin- spin interac- 
tion Hj is present when the ODF laser beams are on. We choose 
jUR — (Oi = nx In/Xarm SO that for small detunings from the COM 
mode (coi), the spins are decoupled from the motion at the end of 
the drive, (b) Evolution of a single spin prior to the application of the 
spin-echo 7i-pulse. c,d) Plots exhibiting spin precession oc (a^) due to 
Hi as a function of initial tipping angle 9i . The error bars are statisti- 
cal (s.d., /I ~ 200). The plots show typical experimental data (black) 
and single-parameter fits to Eq|4] For an antiferromagnetic coupling 
(AFM) and Xarm 250 jus, /jr - (Oi = 2% x 4.0 kHz, Ir = 1 .4 W-cm"^ 
we obtain J/Ir = 27C x 25 Hz • • cm^ (yellow fit). Longer drive 
periods and higher laser intensity Ir yield a larger precession. For 
^arm = 350^s, jiiR - (Oi = 2% X 2.9 kHz, Ir = 1.9W-cm~^ we ob- 
tain ///| = 271 X 55 Hz • • cm^ (red fit). The data in this plot 
is typical of the experiments conducted for benchmarking, (d) For a 
much stronger interaction, Eq.|4]cannot be used to obtain / since the 
MF assumption is no longer valid (see Supplementary Information). 
Also, here we used a small negative detuning (CO2 <C a^/? < coi) which 
gives a long-range ferromagnetic (FM) interaction. For these exper- 
iments we set cOr = 271 X 45.6 kHz. (e) Benchmarking results for an 
ion-crystal with N — 206 ±10 ions. Each point is generated by mea- 
suring / as in c) and measuring the laser beam intensity at the ions 
Ir. The vertical error bars are dominated by uncertainty in Ir (see 
Supplementary Information). The solid line (red) is the prediction 
of MF theory that accounts for couplings to all N transverse modes; 
there are no free parameters. The line's breadth reflects experimental 
uncertainty in the angle 9/^=4. 8°±0. 25. The MF prediction for the 
average power-law scaling coefficient a is drawn in green (right axis, 
linear scale). 
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state, Bj gives rise to spin precession about z in excess of that 
expected due to simple Larmor precession. The experiment 
sequence shown in Fig. measures this excess precession, 
averaged over all spins in the crystal. At the outset, each spin 
is prepared in followed by a rotation about the x-axis by 
an angle 0i . ^/ is applied during the arms of a spin-echo; pre- 
cession (a^) coherently adds over both spin-echo intervals. 
The final 7r/2-pulse maps precession out of the initial jz-plane 
into excursions along z (above or below the equatorial plane 
of the Bloch sphere) that are resolved by projective spin mea- 
surement along t 

As a function of "tipping-angle" 0i we detect global state- 
dependent fluorescence bright, ||) dark) using a photo- 
multiplier tube. This measurement permits a systematic study 
of MF-induced spin precession averaged over all particles, 

^l^pj = 2 (^L^i:^,,/^) COS01 = 2/COS01. The prob- 
ability of detecting |t) at the end of the sequence is 

P(|t)) = ^ (1 + exp(-r- 2x^,^) sin(0i) sin(2/cos(0i) • 2x«,^)) , 

(4) 

and a single-parameter fit to experimental data yields /. De- 
coherence due to spontaneous emission is accounted for by 
r and is fixed by independent measurement of the ODF laser 
beam intensities Ir (F oc see Supplementary Information). 

In Fig. |3]:,d we show representative measurements of ex- 
cess precession due to Hj for different values of spin cou- 
pling strength (determined by Jtj Ir) and interaction dura- 
tion 2Zarm' The cxccss spin precession varies periodically with 
01 (period n) and larger interaction strength yields greater pre- 
cession, manifested in our experiment as a larger amplitude 
modulation of P(|t)). Data agree with Eq.|4]and allow direct 
extraction of / for given experimental conditions. In Fig. [S^ 
we plot /, normalized by /| (Ir independently measured), as 
function of detuning jur — coi (N = 206 ±10 ions). Using no 
free parameters, we find excellent agreement with / obtained 
by averaging over all Jij, where the Jij were calculated by 
including couplings to all N transverse modes (Eq.[3]). 

The MF interpretation of our benchmarking measurement 
tolerates only weak spin- spin correlations. Therefore, in 
the benchmarking regime we apply a weak interaction (/• 
2Xarm <C /4, scc Supplementary Information). In a quan- 
tum simulation, the same interaction is applied at greater 
power producing quantum spin- spin correlations. In the 
present configuration of our apparatus, spontaneous emission 
due to the ODF laser beams (F in Eq.|4]) is the dominant source 
of decoherence. With modest laser powers of 4mW/beam and 
a detuning /i/^ — coi = 271 x 2kHz, we obtain J^2nx 0.5 kHz 
and T/Jr^ 0.06. The expected spin-squeezing {J^) due to this 
interaction is 5 dB, limited by spontaneous emission. The ratio 
r//can be reduced (/increased) by a factor of 50 by increas- 
ing %R to 35°, a likely prerequisite for access to the shortest- 
range, dipole-dipole coupling regime (a ^ 3). At present, ge- 
ometric constraints limit 0/? < 5°; we plan upgrades to our 
apparatus to permit %r = 35°. We also note that relaxation of 
the constraint ~ —F^ can also reduce F// 

In summary, this work establishes the suitability of a Pen- 
ning trap apparatus to pursue quantum simulation in a regime 



inaccessible to classical computation. Our approach is based 
on naturally occurring 2D Coulomb (Wigner) crystals with 
hundreds of ion qubits, a novel experimental system that does 
not require demanding trap-engineering efforts. Experimen- 
tally, we used an optical-dipole force to engineer a tunable- 
range spin- spin interaction and benchmarked its interaction 
strength. Excellent agreement is obtained with the predictions 
of mean field theory and atomistic calculations that predict 
a power-law antiferromegnetic spin coupling Jij oc Jrj^ for 
0.05 <a< 1.4. 

With this work as a foundation, we anticipate a variety of 
future investigations. For example, simultaneous application 
of non-commuting interactions Hb and Hj is expected to give 
rise to quantum phase transitions; Hj may be antiferromag- 
netic (jLiR > COi) or ferromagnetic (CO2 <C /i/? < COi)). Geometric 
modifications to our apparatus will permit access to larger Br 
and antiferromagnetic dipole-dipole-type couplings (a ^ 3). 
Improved image processing software will permit direct mea- 
surement of spin- spin correlation functions using our existing 
single-spin-resolving imaging system (Fig.[T]). 

METHODS 

In a frame rotating at frequency 0)^, the trap potential is 

^(^(r,z) = ^Mco?(z' + Pr2), (5) 

where q is the ion charge, M is the single-ion mass and 
P = co^co^^(^lc — ^r) — 1/2. The ^Be+ cyclotron frequency 
is Q.C = Boq/M = 2% X 1 .6 MHz and the frequency of the 
ions' harmonic center-of-mass (COM) motion along £ is coi = 
27r X 795 kHz. Ion rotation is precisely controlled with an ex- 
ternal rotating quadrupole potential | 21 1. For 100 <N < 300 
we set co^ ~ 271 X 45 kHz so that the radial confinement is suf- 
ficiently weak that a cloud of ions relaxes into a single 2D 
plane (P <C 1). Upon Doppler laser cooling the ions' motional 
degrees of freedom (Tcom ^ 1 mK) 1301 . the ions naturally 
form a 2D Coulomb crystal on a triangular lattice, the geom- 
etry that minimizes the energy of their mutual Coulomb po- 
tential energy. The crystal has N transverse eigenmodes co^ 
with eigenfunctions 5^; the COM mode COi is the highest fre- 
quency mode (see Fig.|2]). 

The spin-dependent ODF is generated by a pair of off- 
resonance laser beams with angular separation Qr ~ 4.8° and 
difference frequency jur (see Fig.[T]). The result is a traveling 
ID optical lattice of wavelength .7 ^m whose 

wavefronts propagate along £, traversing the ion crystal at fre- 
quency iHR/2n. Alignment of is crucial for proper spin- 
spin coupling (see Supplementary Information). The lattice's 
polarization gradient induces a differential AC Stark shift on 
the qubit states (a spin-dependent force). We choose operat- 
ing conditions that give F^ ^ —F^, where F^ = Focos{iLiRt)t 
For reference, if the single-beam intensity at the ions is Ir = 
1 W • cm-^, we obtain - 1.4 x 10"^^ N. 

Small coherent displacements that produce negligible spin- 
motion entanglement (as required by Eq. |3]) are obtained for 



detunings satisfying 

fi \iAR - co^ I > Fq ^fi{2nm + 1 ) / (2Mco^) . (6) 

This is a more stringent criterion than that used by others 
|[T3l |29l as it includes an additional to account for a 
typical distribution of composite spin states. Moreover, we 
also include a correction factor for finite temperature ^ 

kBT/h(Om. 
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SUPPLEMENTARY INFORMATION 

A review of ion confinement in Penning traps and discus- 
sion of a variety of equifibrium states, including 2D Coulomb 
crystals (ionic Wigner crystals), can be found in 1 211 [3111321 . 
A number of authors have theoretically analyzed and dis- 
cussed the prospects of using 2D Coulomb crystals for quan- 
tum information and computation (281 [33^37 J . The engi- 
neered Ising interaction, which we report here, builds on our 
previous experimental work using 2D Coulomb crystals for 
high-fidelity quantum control I26,.38.-.41 1. Below we discuss 
some details of this new capability. 



I. SPIN INITIALIZATION, CONTROL, AND 
MEASUREMENT 

Reference f26i\ gives a detailed discussion of our spin ini- 
tialization, control, and measurement capabilities with planar 
ion arrays in Penning traps. Here we briefly summarize some 
of that discussion, emphasizing aspects relevant for the mea- 
surements reported here. Figure |4] shows the relevant ^Be+ 
energy levels. We use the valence electron spin states paral- 
lel It) = |m/ = + and antiparallel ||) = |m/ = — ^) to the 
applied magnetic field of the Penning trap as the two-level 
system or qubit. In the 4.46 T magnetic field of the trap, 
these levels are split by approximately = 2nx 124 GHz. 
The ^Be^ nucleus has spin 7 = 3/2. However, we optically 
pump the nuclear spin to the m/ = +3/2 level |42 |, where it 
remains throughout the duration of an experiment. The ions 
are Doppler laser-cooled to a temperature ~ 1 mK 1 30] by a 
1 MHz line width, 313 nm laser tuned approximately 10 MHz 
below the |t) p^3/2^/ = +3/2) cycling transition. Spins 
in the ||) state are efficiently optically pumped to the |t) state 
by a laser tuned to the ||) p^3/2^/ = +1/^) transition. 
The repump beam and the main Doppler laser cooling beam 
are directed along the magnetic field (£-axis). Powers are a 
few milli-Watt with laser beam waists of ~ 1 mm. In ad- 
dition, a weak Doppler laser cooling beam ('^40 ^m waist) 
directed perpendicularly to the £-axis directly Doppler cools 
the perpendicular degrees of freedom. A typical experimental 
cycle starts with 10 ms to 20 ms of combined Doppler laser 
cooling and repumping. The repump laser remains on for an- 
other 3 ms after the Doppler cooling laser is turned off. The 
fidelity of the |t) state preparation is estimated to be very high 
(> 99.9%) Ei. 

Low-phase-noise microwave radiation at 124 GHz is used 
to rotate the spins through the magnetic dipole interaction 

Hb = gjuslli^jjit) • , where ^^(0 is the applied mi- 

crowave field (predominantly perpendicular to z), g — 2 is the 
electron g-factor, and jub is the Bohr magnetron. The fidelity 
of a 7i-pulse was measured to be greater than 99.9% in a ran- 
dom benchmarking experiment 1 26 1 . The microwave source 
consisted of an agile 15.5 GHz source followed by an ampli- 
fier multiplier chain with 150 mW output power at 124 GHz. 
The 15.5 GHz source is obtained by mixing, with a single- 
side-band mixer, the output of a 15.2 GHz dielectric resonator 



oscillator (DRO) with the output of a 300 MHz direct digi- 
tal synthesizer (DDS) that is under field programmable gate 
array (FPGA) control. On/off switching of the 124 GHz mi- 
crowaves is done at 15.5 GHz before the amplifier multiplier 
chain. The microwaves are transported to the ions down the 
bore of the magnet with a rigid waveguide and directed onto 
the ions with a horn located between the ring and endcap elec- 
trodes of the trap. With this arrangement the microwave hard- 
ware does not block optical access along the magnetic field 
axis, enabling imaging of the ion resonance fluorescence scat- 
tered along the magnetic field (top-view image — see [Hi]). 
We obtain TT-pulses of 70^s duration with the 150 mW out- 
put power of the amplifier multiplier chain and the setup de- 
scribed here. The measured spin-echo coherence duration (T2) 
is ~ 100 ms. 

At the end of an experimental sequence we turn on the 
Doppler cooling laser and make a projective measurement of 
the ion spin state through state-dependent resonance fluores- 
cence. With the Doppler cooling laser on, an ion in the |t) 
state scatters ~ 10^ photons/s while an ion in ||) is dark. For 
the spin precession measurements reported here we performed 
a global spin-state detection. Specifically we detected, with 
f/5 light collection and a photomultiplier tube, the resonance 
fluorescence from all the ions in a direction perpendicular to 
the magnetic field (the side-view). For detection periods of 
~ 50 ms the detection fidelity is high, typically limited by 
quantum projection noise. Here we used short detection pe- 
riods of ^ 500/is, from which we would detect ~ 1 photon 
for each bright state Typically, each experiment was re- 
peated ~ 100 times and averaged, resulting in a few percent 
uncertainty due to shot noise in the measurement of P{'t). 

The spin-precession signal used to benchmark spin-spin 
coupling in the manuscript relied on a global spin-state mea- 
surement via side-view fluorescence collected on a photo- 
multiplier tube. In the future, we anticipate that time-resolved 
top- view images such as that shown in Fig. 1 will be used to 
obtain the spin state of individual ions. As discussed in the 
manuscript, ion crystal rotation at co^ is phase-locked locked 
to an external oscillator. We use an imaging photomultiplier 
tube to record (x^y^t) for each photon. Rotating-frame im- 
ages are generated computationally given (x^y^t) and co^, a 
technique established in 2001 1431 . In linear Paul trap ex- 
periments, determination of ion spin-state is possible with as 
few as 10 photons/ion [441 . At present we await the arrival 
of a new (x^y^t) detector system capable of a detection rate of 
5 X 10^ Hz. We anticipate this will enable high-fidelity spin- 
state measurement of a 300-ion crystal in 10 ms. From a 
suite of identical experiments, the spin- spin correlation func- 
tion can be computed. We believe that our ability to resolve 
single ions, even in the presence of rotation at CO/? indicates 
a path forward in performing individually resolved measure- 
ments of fluorescence correlations between ions. 

Previous measurements have elucidated some of the possi- 
ble limiting mechanisms. For instance, the stability of crys- 
tal orientation in the rotating frame was studied for spheri- 
cal crystals of ~ 10,000 ions. In these experiments the ori- 
entation was observed to precess uniformly for durations of 
> 10 s, then suddenly slip by a large angle before again re- 
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suming slow precession f4T|. This "stick-slip" motion can 
be followed and easily corrected. Figure la of the manuscript 
was generated from 600 s of integration after stick- slip correc- 
tions. Other potential issues include ion loss and background 
gas collisions. As the trap depth is 1 eV, background gas 
collisions do not result ion loss. Collisions with background 
hydrogen generate BeH+ (1 per 6 minutes for A/^ ~ 300 ^Be+ 
ions) which collect at the crystal perimeter due to centrifu- 
gal separation. These effects do not appear to limit top- view 
imaging fidelity, as demonstrated in Fig. la of the manuscript. 
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II. OPTICAL-DIPOLE-FORCE LASER SETTINGS 

Figure |5] shows a simple sketch of the optical dipole force 
(ODF) laser beam set-up. As discussed below, the frequency 
as well as the beam polarizations were chosen to null the 
AC Stark shift from an individual beam and to produce a state- 
dependent force that is equal in magnitude but opposite in sign 
for the It) and ||) qubit states (F^ = —F^). This reduced the 
system sensitivity to laser intensity fluctuations. For exam- 
ple, if F^ / —F^, then the interaction induced by the optical 
dipole force will include terms linear in the dj's. These terms 
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Figure 4. Relevant energy levels of ^Be+ at Bq = 4.46 T (not drawn 
to scale). We show only m/ = + ^ levels that are prepared experimen- 
tally through optical pumping. The ^»S'i/2 ~ ^^3/2 transition wave- 
length is ~ 313 nm. A resonant laser beam provides Doppler laser 
cooling and state discrimination; a second laser beam repumps ||) 
to It). The optical dipole force (ODF) interaction is due to a pair 
of beams (derived from the same laser) with relative detuning /jr. 
The qubit splitting AEq/H ^2nx 124 GHz. A low-phase-noise mi- 
crowave source at 124 GHz provides full global control over spins. 



Figure 5. Sketch of optical dipole force (ODF) laser beam setup, 
(a) The ODF laser beams lie in the y-z plane at angles ±6/?/2 with 
respect to the j-axis. (b) View looking in the —y direction. The 
beams are linearly polarized but with different polarization angles 
relative to vertical polarization. 



can be canceled with spin-echo techniques, but this requires 
that laser intensity fluctuations are small. Likewise, by ad- 
justing the laser polarization to null the AC Stark shift from 
a single beam, we mitigated qubit decoherence due to laser 
intensity fluctuations. For the benchmarking measurements 
described here we did not actively stabilize (i.e., noise eat) the 
laser beam intensity. 

The off-resonant laser beam frequency was detuned from 
the cycling transition (|t) ^ IP3/2, = 3/2)) by Ar c:^ 
—63.8 GHz. This gives detunings of +15.6 GHz and 
-26.1 GHz, respectively, from the |t) \p^i2,mj = 1/2) 



and II) 1^3/2, = —1/^) transitions (Fig.|4|. Laser beam 
waists were w^ — l lO/^m in the vertical (z-direction) and ^ 
1 mm in the horizontal direction. Here we define the waist as 
the distance from the center of the beam over which the elec- 



tric field intensity decreases by l/e^ (i.e. I{z) ~ e 



). 



With the small 2.4° incident angle each beam makes with re- 
spect to the plane of the crystal, this provided a uniform elec- 
tric field with < 10% intensity variation across ion crystal ar- 
rays with N < 250. 

The ODF laser beams were linearly polarized at nonzero 
angles with respect to the z-axis. Let 



Eu (r, t) = £[/£'[/ cos {^u ' ^ — ^ut^ 



(7) 



denote the electric fields of the upper and lower ODF beams. 
If (\>p is the angle of the laser beam electric-field polarization 
with respect to vertical polarization (e-x = 0), then the AC 
Stark shift of the qubit states when illuminated by a single 
beam can be written 



^t.acss = A|COS^ {(^p) +5|Sin^ {(^p) 



(8) 



where A|(A;) is the Stark shift of the |t)(||)) state for a %- 
polarized beam (8 parallel to the z-axis) and B^(B^) is the Stark 
shift of the |t)(|i)) state for a a-polarized beam (8 perpendic- 
ular to the z-axis). (Here we neglect the small a polarization 
(oc sin 2.4'') that exists when ^p = 0.) The Stark shift of the 
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qubit transition is 

^acss = (At - A;) cos^ ((^^) + {B^ - B^) sin^ {(^p) . (9) 

If —A^ and B^ — B^ have opposite signs, there is an angle 
which makes A^^^^ = 0. For a laser detuning of Ar = —63.8 
GHz, A^,,, = 0at^^^±65^. 

With Aacss = for each ODF laser beam, we exploit the 
freedom to choose their polarization in order to obtain a state- 
dependent force. Specifically, we choose Eu to have a polar- 
ization given by (^p^u = 65^ and to have a polarization given 
by = —65^. In this case the interference term in the ex- 
pression for the electric field intensity (^Eu + E^^ produces a 

polarization gradient which results in spatially dependent AC 
Stark shifts 



(A^cos^ {^p) —B^ sin^ i^p)) 2cos(6^- 
(A; cos^ {(^p) - B^ sin^ {(^p)) 2 cos (5^ • 



(10) 



ku-ki 



2/: sin 



in(^)i 



is the 



for the qubit levels. Here 5^ ; 

wave vector difference between the two ODF laser beams, 
j^R = — is the ODF beat note, and (^p = \(^p^u\ = \^p,i \ • 
The spatially dependent AC Stark shift produces a state- 
dependent force = ^ot,^ (^^ ' ^ ~ I^R^) where 



= -2dk (A^cos^ {(^p) -B^sin^ {d^p)) 
Foi = -28k (A;cos2 {(^p) -B^sin^ ((^^)) . 



(11) 



In general F^^ 7^ —FqI- We operate at A/^ = —63.8 GHz where 
for Aacss = we also obtain F^^ = —F^i = Fq 

For a given ^p^u , ^p,i , and Ar we use straightforward atomic 
physics along with well known values for the energy levels 
and matrix elements of ^Be+ to calculate Fo as a function of 
the electric field intensity = ^ = ^ |^[/ P at the cen- 
ter of the laser beams. For %r = 4.8° and Ir = I W/cm^ , 
Fo = lAxlO-^^ N. 

Stronger forces can be generated after experimental mod- 
ification to our apparatus to permit Qr ^ 35^ . At this an- 
gle larger detunings jur — COi are required to satisfy Eq. 6 in 
Methods. With our definition Hj = ^Y.i<jJij^]^)^ the in- 
teraction strength between two ions / and j is Jij/N. Con- 
sider for example N — 211 ions, co^ = 2:r x 45.6 kHz, pow- 
ers of Ir = 20 mW/beam (12.5 W/cm^) and Qr = 35^ The 
spin-motion entanglement constraint is satisfied by a detun- 
ing of /i/? — COi = 271 X 100 kHz. In this case we obtain 
Jij/{2KN) - {560 H z) {do /dijY-\ where do - 20iLim is the 
typical nearest neighbor separation. 



III. WAVEFRONT ALIGNMENT 

The ODF laser beams produce a one dimensional (ID) op- 
tical lattice characterized by the effective wave vector dk and 
beat note jur. In Sec.|ll]we assumed that dk\\t or equivalently 
that the wavefronts of the lattice were aligned perpendicular 
to the z-axis (magnetic field axis). If the wavefronts are not 
normal to the z-axis as sketched in Fig. [6j then the time de- 
pendence of the ODF seen by an ion in the rotating frame 




do-400 |um ' 



Figure 6. (a) Sketch of the ODF laser beam geometry used to gen- 
erate the one dimensional (ID) traveling optical lattice, (b) Sketch 
of the ID optical lattice wave fronts (red lines). These wave fronts 
need to be aligned with with the ion planar array (represented by the 
blue dots). Here Xr = 27i/|A^| ^ 3.7 /jm and Qerr denotes the angle 
of misalignment. 400 ^m is the typical array diameter for N ^ 200 
ions. With the wavefront alignment technique discussed in the text 
we obtain Oerr < 0.05°. 



depends on the {x^y) position of the ion. This complicates the 
effective spin- spin interactions generated by the ODF but can 
be adequately mitigated by careful alignment. 

Alignment of the ODF laser beams is obtained by a tech- 
nique that makes top-view images (images of the ion reso- 
nance fluorescence scattered along the magnetic field) from 
a single plane of ions sensitive to ODF wave front misalign- 
ment. For this measurement we set jlir = (stationary ID lat- 
tice) and detune the frequency of the ODF laser beams approx- 
imately 0.5 GHz below the |t) ^ ^^3/2 m/ = +3/2) Doppler 
cooling transition. This small detuning generates sufficiently 




Figure 7. Top- view image of the spatially inhomogeneous fluores- 
cence from a single ion plane produced by the AC Stark shift from a 
static {/JR = 0) ODF lattice with misaligned wave fronts. Dark bands 
are regions of high standing wave electric field intensity (parallel 
to the dashed yellow line). The bright horizontal feature bisecting 
the center of the image is fluorescence from the weak Doppler laser 
cooling beam directed perpendicular to the magnetic field. The im- 
age was obtained by subtracting a background image with the ODF 
beams off (1 s integration). 
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obtained by expanding Hi to first order in 8a^- = d^j — / d^- 



Figure 8. (a) Sequence used to measure spin precession due to the 
effective mean field generated by the engineered Ising interaction. 
The spins are initially prepared in |t). The first pulse rotates the 
spins by an angle 9i about an axis (defined to be the i-axis in the 
rotating frame) in the equatorial plane. The subsequent 7i and 7i/2 
pulses occur about the j-axis. (b) Evolution of a single spin prior to 
the application of the spin-echo 7i-pulse. 



large AC Stark shifts on the cooling transition to measurably 
change the ion scatter rate from the Doppler cooling laser. 
With the Doppler cooling laser turned on and the ODF beams 
turned off, we observe a spatially uniform, time-averaged im- 
age of a rotating planar crystal. With the ODF beams on, ions 
located in regions of high electric field intensity at the anti- 
nodes of the optical lattice are Stark- shifted out of resonance 
with the Doppler cooling laser. This is the cause of the dark 
bands in the top- view image shown in Fig. [7] We adjust the 
ODF laser beams based on this real-time imaging to optimize 
their alignment. Improved alignment is indicated by a fringe 
pattern of longer wavelength. With this technique we have 
aligned the ODF wave fronts with the planar array to better 
than -0.05°. 

The image in Fig. [7] is typical of what we obtain with 1 s 
integration duration. This indicates the ID lattice was stable 
during the integration period and shows the phase stability of 
our ID lattice of better than 1 s. 

We note that direct fluorescence imaging of the ID lat- 
tice, for example by tuning the ODF laser resonant with the 
Doppler cooling transition, is not viable. Even at low powers, 
resonantly scattered photons across the large horizontal waist 
of the ODF beams apply a large torque, causing the rotation 
frequency and radius of the array to rapidly change, typically 
driving the ions into orbits of very large radius. 

We have also used phase coherent Doppler velocimetry to 
improve the ODF wavefront alignment ll45ll . But the top-view 
imaging technique discussed here and shown in Fig. |7] pro- 
vides more information on the angle and direction of mis- 
alignment, which greatly improves the ODF-crystal alignment 
process. 



IV. MODELING MEAN FIELD SPIN PRECESSION 



For any spin j. 



con- 



0, which implies . 

stant under application of Hj. The Hamiltonian Hf^^ describ- 
ing the mean field response of spin j to the application Hi is 



L<,.y,-j((af)5^;. + (d})5^;) 



Defining Bj = ^ ^i^jJij (af ), Hf^ can be written 



(12) 



(13) 



The mean field Heisenberg equations of motion for spin 
is obtained from the commutator of and Hf^^ . However, 

because 6a ^ = al — (d^) and (d^) is a constant, it is clear 
that Eq. [T3[describes spin precession about the z-axis at fre- 
quency B^. Our observable is the spin precession (in excess 
of ordinary Larmor precession) averaged over all the spins 

We use the spin-echo sequence in Fig. [8] to measure a pre- 
cession proportional to the expectation value of the spin pro- 
jection along z ((6^)). The spin-echo sequence minimizes 
contributions to spin precession that are not ^^(g^). Specifi- 
cally, the spin echo cancels a constant spin precession inde- 
pendent of (d^) (e.g., due to slow uncontrolled magnetic field 
fluctuations), but precession o^(a^) coherently adds in the two 
arms of the sequence. 

The first pulse sets (dj) = cos 0i . The interaction Hj is then 
applied by turning on the ODF laser beams for a period Xarm- 
During this interval (cj) is a constant and mean field theory 
predicts an average spin precession angle of 2/cos(0i) • Xarm- 
The TT-pulse changes (dj) — cos(0i) and the spin preces- 
sion angle 2/cos(0i) -Xarm to — 2/cos(0i) -Xarm- At the end 
of the second Hj interaction of duration Xarm, the total preces- 
sion angle of the spins is — 2/cos(0i) • 2Xarm- The final 7z/2- 
pulse is about an axis shifted by 90^ from the first 0i -pulse. 
This pulse converts precession out of the initial yz plane into 
excursions above or below the equatorial plane of the Bloch 
sphere, which we measure. 

The evolution operator of the measurement sequence Useq 
is obtained in a straight forward manner from the individual 
evolution operators from each segment of the sequence. 

Use, =^{y^l)- UiHD • R (j, n) • U {H^) • ^ (x, 01 ) . 

(14) 

cos(0i/2) -/sin(0i/2) 
-/sin(0i/2) cos(0i/2) 
1 -1 
1 1 



Here ^(x,0i) 



-1 

1 



and^(j,|) 



2 



lution is given by 



. The mean field evo- 



u{hD = 



QXp{-iJ{Cz)%rm) _ 

QXp{iJ{az)Xarm) 



(15) 



where (d^) = cos(0i) in the first arm of the spin-echo se- 
quence and (d^) = — cos(0i) in the second arm. At the end 
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Figure 9. We determine from the AC Stark shift of the qubit tran- 
sition as a function of ODF laser beam polarization angle The 
triangles and filled circles are separate measurements done on the 
upper and lower ODF beams. 1r is determined independently for 
each ODF beam by a fit to Eq.|9] The error bars are smaller than the 
plot points. 



of the sequence we detect the |t) state probability. Explicit 
computation gives 



Kti^/..,it>|' 



1 



{ 1 + sin (01 ) • sin [2/cos (0i ) • lXarm\ } 



(16) 



V. OPTICAL DIPOLE FORCE LASER INTENSITY 
CALIBRATION 

The spin-spin couplings (/) measured here depend on the 
square of the ODF electric field intensity Ir. Therefore, care- 
ful calibration of the electric field intensity was important in 
benchmarking the strength of the interactions with mean-field 
predictions. We separately determined Ir for each beam by 
measuring the AC Stark shift of the qubit transition as a func- 
tion of the laser beam polarization angle (|)p, as shown in Fig. 
[9] The qubit transition frequency was measured by fitting for 
the center frequency of a Rabi-resonance profile. The po- 
larization angle was varied by rotating a \ll plate. The 
Stark shift measurements were fit to Eq. |9] which provided 
values for A and B (Stark shifts for TT-polarization and a- 
polarization). Values for A and B were then used with atomic 
physics calculations to determine Ir, the electric field intensity 
of the ODF laser beams, with a fractional uncertainty of ~ 5%. 
Frequent intensity calibration measurements were taken dur- 
ing a benchmarking run. Slow drifts to the laser intensity be- 
tween calibrations added another 5% uncertainty, which we 
add in quadrature to the fitted uncertainty. 



VI. SPONTANEOUS EMISSION 

Decoherence due to spontaneous emission has been well 
studied in this system [41 1. The qubit levels are closed under 
spontaneous light scattering; that is, spontaneous light scat- 
tering does not optically pump the ion to a different ground 



state level outside of the two qubit levels. We measure 
spin precession to benchmark the Ising interaction couplings. 
Equivalently we measure the evolution of off-diagonal coher- 
ences between the \X) and ||) levels. Spontaneous-emission- 
induced decay of these coherences is accurately modeled by 
Eq. (8) of Ref. 1411 . We add this time dependence to the uni- 

This results in 



IV 



tary evolution JJ {flf^^^ discussed in Sec. 
the expression 

P (It)) = ^ (1 + exp (-r • IXarm) sin (01 ) sin (2/cos (0i ) • IXarm)) , 

(17) 

where P{W)) is the probability of detecting a spin in W) 
at the end of the spin-echo sequence. The only difference 
with Eq.[T6]is the factor exp(— F • 2Xarm) due to the decay of 
the off-diagonal elements of the density matrix. Here F ac- 
counts for decoherence due to spontaneous emission. From 
im, F = ^ iXRam^^el) h^s Contributions from both Raman 
scattering and elastic Rayleigh scattering. A straightforward 
atomic physics calculation relates F to the ODF laser beam 



intensity Ir. In fits of Eq. 17 to the spin-precession measure- 
ments, we fix F at the value determined by the laser intensity 
calibrations. For ^R = -63.8 GHz, (^p = ±65.3^, and Ir = 1 
W/cm^, we calculate F = 82 s~^. 

With our present set-up we can achieve /::^ F for detunings 
\liiR — (^i \ ^ 10 kHz. For these detunings it will be possible to 
simulate quantum effects beyond mean field theory such as 
spin squeezing and spin-depolarization due to many-body in- 
teractions. For example we calculate that with 4 mW/beam 
and an ODF beatnote detuning of |^/? — COi | =2 kHz , / ^ 
271 X 0.5 kHz and F//~ 0.06. The potential spin squeezing 
due to this interaction is 5 dB, limited by spontaneous emis- 
sion. For many simulations it will be desirable to achieve 
J/N r where / is a nearest-neighbor coupling strength. 
(J ^ J for small detunings \jUR — coi |.) The most straight for- 
ward strategy to achieve this condition in our set-up appears 
to be to increase the angle Qr between the Raman beams. The 
ratio //F scales as sin^ (0/^/2). Therefore decoherence due to 
spontaneous emission can be dramatically reduced with an or- 
der of magnitude increase in 0/^. We note that a large increase 
in dR will require implementation of a sub-Doppler cooling 



scheme to remain in the Lamb-Dicke regime (see Sec. VII). 

Different ODF laser detunings can likely help reduce the 
impact of spontaneous emission (with the complication that 
^ —F^), but very large laser detunings A obtained by tun- 
ing the ODF laser beam frequencies outside the Pi /2 — ^3 /2 
manifold appear unlikely to help. This is because both the 
interaction strength and spontaneous emission scale as 1/A^ 
for our qubit, which is not a clock transition. For trapped-ion 
experiments in low magnetic field, the impact of spontaneous 
emission is typically minimized by tuning the ODF laser beam 
frequencies between the P1/2 and manifolds. Because the 
^Be+ P- state fine structure (197 GHz) is comparable to the 
Zeeman splittings for our magnetic field (4.5 T), this approach 
does not significantly help. However, it can help with heavier 
mass ions (e.g., Mg+). Reference (TP\ discusses a number of 
strategies that can minimize the impact of spontaneous emis- 
sion when the P-state fine structure is large compared to the 
Zeeman interaction energy. 
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VII. LAMB-DICKE CONFINEMENT 

Throughout this manuscript we impHcitly assume 
that the appHed state-dependent force, = 
Fo^^lcos{dk- z — jLiRt), is constant across the spatial ex- 
tent of an ion's wave function. Specifically, for a planar array 
located at z = 0, we assume the force on any ion is given by 

F^^^{z = 0,t)=Foucos{iLiRt) . (18) 

The extent to which this is true is quantified by the parameter 

r{ind,i = ^k-Zrms4 (19) 

where Zrms.i is the root mean square (rms) axial extent of the 
wave function of ion / and 6k = 2n/XR, where Xr = 3.7 jum for 
Qr = 4.8^. We note that r\ind,i is a Lamb-Dicke confinement 
parameter for a individual ion (ind), not to be confused with 
the usual Lamb-Dicke parameter that is defined in terms of 
the ground-state wave function. In the limit that the Coulomb 
interaction energy between ions is a small perturbation to the 
axial potential of the trap, we can think of each ion as a single 
ion confined in the external trap potential, and approximate 

Zrms.i — ^ 2Mo^ ' 1, where M is the mass of an indi- 

vidual ion. Recently we have completed careful axial tem- 
perature measurements for single-plane ion crystals 1301 . We 
find ~ 1 mK for the axial COM mode and ~ 0.4 mK for 
the higher order transverse modes. For co^ = 271 • 795 kHz 
and conservatively assuming T = 1 mK, we calculate n ~ 26, 
Zrms.i - 190 nm, and \\ind,i ^ 0.32. 

For the benchmarking measurements described here we 
typically set the rotation frequency of the ion array about 
^0.5 kHz below the rotation frequency of the 1 ^ 2 plane 
transition, determined experimentally from side-view images 
f21 1. In this case, the spectrum of the transverse (axial) modes 
is broad, and we underestimate Zrms.i in the above analysis. An 
improved estimate of Zrms.i is obtained by summing the con- 
tributions from all of the transverse modes m 



Zrms,i — f T^, ipi.m) 



fl 



■(2^^ + l) 



1/2 



(20) 



where Rm — ksT /h(})ni is the mean thermal occupation of 
mode m. We assume here that every mode is characterized by 
the same temperature T. For N = 211 and co^ = 271 x 795 kHz, 
the rotation frequency of the 1 ^ 2 plane transition is co^ = 
27r X 46.1 kHz. For CO/^ = 27r x 45.6 kHz, the spectrum of the 
axial modes ranges from 795 kHz down to 224 kHz. With 
Eq.[20|we calculate Zrms.i — 520 nm in the center of the array, 
decreasing to Zrms.i — 250 nm at the array edge, corresponding 
to ^ind.i — 0-89 for / near the array center and r\ind,i = 0-42 for 
/ near the array boundary. 

The wavefront alignment discussed in a previous section 
can be viewed as achieving a type of Lamb-Dicke confine- 
ment. Let derr dcnotc the angle of misalignment between the 
planar array and the ODF ID lattice (Fig. [6]). Relative to the 
ODF lattice wavefronts, the rotation of the array produces a 
time-dependent shift in the axial position of an ion that can 



be written as z{t) = /?sin(0g^^) sin(co^r + (p). Here R is the 
distance of the ion from the center of the array, and (p is deter- 
mined by the azimuthal position of the ion in the array. The 
force on this ion is then 

F^^^{t)=Foucos{dkRsmeerr-^^ri{o}rt^(?)-ILiRt) . (21) 



For Eq. [2T]to approximate Eq. 1 8 we desire dk 2Rp sin derr < 
1, where Rp is the array radius. For Rp 200 /^m (typical for 
N = 200) and derr — 0.05° (see Wavefront alignment section), 
we calculate 5^ • 2Rp sinO^^^ ^ 0.6. 



VIII. ION LATTICE CONFIGURATION AT 
EQUILIBRIUM AND TRANSVERSE NORMAL MODES 

The ion equilibrium positions fi at zero temperature are cal- 
culated by minimizing the Euler-Lagrange action for the Pen- 
ning trap potentials and the constraint that the ions lie in a 
plane (at z = 0). The solution is a triangular lattice with a 
lattice constant that increases as one moves radially outward, 
that has a smooth unfaceted edge and exhibits a degradation in 
orientational order near the crystal perimeter. The transverse 
(along z) phonon modes (o)^, bm) are obtained by Taylor ex- 
pansion of the potential about the ion equilibrium positions 
In-plane modes (along x^y) can also be calculated, by solv- 
ing a quadratic eigenvalue problem (due to inclusion of the 
centrifugal and Coriolis forces). The details of the transverse 
mode calculation is discussed in this section. The problem has 
also been solved for longitudinal and transverse modes in ID 

(e.g., mmy 

In general, the Lagrangian for a collection of N ions with 
charge q and mass M in an electromagnetic field —A{fi) is 



L = T-V 



: f I ^Mff - q ((^,- -A{n) • y?j I , (22) 



where = (r/, 9/,^;) is the coordinate of ion /. In a Penning 
trap, the field consists of a uniform magnetic field in the 
^-direction (B = Boz) and a harmonic trapping (anti-trapping) 
electric potential in the z-direction (r-direction) with 
frequency COi . Lastly, an additional external time-dependent 
electric quadrupole potential (the "rotating wall", amplitude 
Vwall at the ions) is applied to control the ion's rotation 
frequency in the trap |48 1. Thus, the scalar potential for ion / 
is 



Q^trap.i + Q^wallJ + Q^Coulomb.i 



= ^Mcof (zf - rf /2) +y^^//rf COS 2(e,- + C0,0 + ^kl^^^kq' 

where the ion-ion separation is rj^k = 0' ~ k=l /47l8o. 

Since B is uniform, the vector potential energy is 



(24) 



The position of the ions at equilibrium is calculated as fol- 
lows. We move to a rotating frame where the ions' coordi- 
nates are stationary by using the coordinate transformation 
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r- = (r-,0-,z-) = (r;,0 + co^r,Zi), a counterclockwise rotation. 
The ions' equilibrium positions r- can be found by solving the 
transformed Euler-Lagrange equations L' . To seed the numer- 
ical solution, we supply an initial guess for the 2D crystal: a 
regular, triangular lattice. Consistent convergence requires ion 
numbers N corresponding to closed shells (e.g., N =121 has 
six closed shells). We find ion equilibrium positions that devi- 
ate from a perfect triangular lattice near the crystal periphery 
and which have an overall ellipticity due to Vwalh 

Given r-, the crystal's transverse eigenmodes can be calcu- 
lated by Taylor expansion of the potential about the equilib- 
rium positions r-. The Lagrangian is 



(25) 



' i=i 



where z- is the axial displacement of /-th ion and Kij = Kji is 
the symmetric stiffness matrix evaluated for the equilibrium 
configuration as 



kct 



n=\ Km I 



(26) 



By minimizing the action 8 J dtL' = with respect to the ro- 
tating frame axial coordinates z-, we obtain N equations of 
motion 



7=1 



■ 0, for / : 



1,2,. 



(27) 



Following standard normal mode analysis, the solution is ob- 
tained by calculating the eigenmodes of the matrix Kij/M. 
The result is eigenvalues co^, for m = 1,2, (with cor- 
responding frequencies (O). 



IX. LIMITS TO THE VALIDITY OF MEAN FIELD 
THEORY 

Here we calculate the domain over which a mean field the- 
ory (MF) treatment of is valid. In the special case of uni- 
form Ising coupling (e.g., < — coi <C coi — CO2), we have 



Hi 




N 



(28) 



where we have neglected a constant offset and = Y4L1 /2 
is the z-component of the total composite spin of the system. 
The interaction strength % is independent of ion-ion separation 
and is given by 



n2Mjji-(oj 



(29) 



We start each experiment with spins all prepared in the 
It) state. Here we use the composite spin picture where this 
state is labeled \J = N/2, Mj = N/2). (Formally \J,Mj) la- 
bels a state that is an eigenstate of P and with eigenvalues 

/(/+ 1) and Mj, where /= Y.i^i/'^ is the total system spin.) 
The first pulse of the spin-precession measurement sequence 
rotates the composite Bloch vector by an angle 0i. After the 
rotation by 9i the state can be written as 



N/2 

|^)= £ C{J,Mj)\J,Mj), 

Mj=-N/2 



(30) 



where the coefficients C{J,Mj) are significantly non-zero for 
a small range AMj centered on Mf\ where ^ ^ cos 61. 
For a coherent spin state, AM j < V^. 

We want to establish that a interaction looks like 

precession oc at least over short durations. A precession 
by an angle ([) about the z-axis is obtained with the operation 

= ^c(/, Mj) e-^^'^ |/, Mj) . (31) 

Mj 

We measure expectation values of the total system spin 7. It 
is sufficient therefore to consider how these expectation 
values transform under the rotation e'^-'^, 

(MfyW) . (32) 

For simphcity, only matrix elements of the total spin in the 
j-direction Jy are considered. Identical expressions are also 
obtained for^. 



= IIm'^^Mj CiJ^ MjfC{J, Mj) {^J, Mj\jy\ J, Mj 
= Lmj {C{J, Mj - l)*C{J, Mj) {J, Mj-l\fy \J, Mj) 

+C{J, Mj + 1 )*C{J, Mj) {J, Mj + l\Jy\J,Mj)}. 

(33) 

The double sum is eliminated by making use of the property 
that Jy has only non-zero matrix elements between states 
with My - Mj = ± 1 . C{J, My ± 1 ) can be defined to be if 
Mj - 1 = -N/2 - 1 or Mj + 1 = N/2 + 1. Similarly 



= ^(•^' ^j)*C{J, Mj)e'«/*e-'^^* (/, M^ |/; | /, Mj 

= li {C{J, Mj - l)*C{J, Mj)e-'^ {J, Mj-\\Jy |7, Mj) 
+C{J, Mj + 1)*C(7, Mj)e'^ {J, Mj + l\fy \J, Mj) } . 

(34) 

This is to be compared with the same matrix element under 
the Hi = interaction, 

{\|/| e'v-'l'Jye~'v-'l* 

> 2 



: Zm'j^^^ CiJ, MjfdJ, Mj)e'^ ' (j, Mj \jy I y, Mj 

= Lmj {c{J, Mj - 1)*C(7, Mj)e-'T(2M,-i)f i^j^ Mj-\\Jy \J, Mj) 

+C(J, Mj + l)*C(y, Mj)e'# (2M/+1> (y, My + 1 1 y; \J, Mj) } 



(35) 
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Equation [35] approximates a rotation about the z-axis (that is, 
approximates Eq. [34]) if 



2X 



t(2Mj-\)^'^t{2Mj^\) 



2X 



(0) 



for all Mj for 



which C{J^Mj) is significantly non-zero. This is satisfied for 
short periods t satisfying '^t -2 (AM /) <C 1. For the initial 
coherent spin state, AM j < \/N, which puts a limit 



5 Particles 




s/N 



(36) 



a) 



on the period t for which the interaction can be 
approximated as a precession oc . The precession 
frequency predicted by the analysis, 

^ • 2Mf'^ 7^ ' ^ — 2% cos 0, is in agreement with the 

precession frequency predicted by the MF analysis. 



Equation[36]provides a limit on the period for which the MF 
analysis is valid. We confirm this limit through explicit cal- 
culation of the spin precession benchmarking sequence for the 
Hi=^J^ interaction with N = 5,50 and 100 spins. The exact 
calculation is possible because the system remains in the sym- 
metric subspace (dimension of Hilbert space is 2N -\- 1). Fig- 
ure [TOja) shows the results of the calculation for five spins. 
Reasonable agreement between the MF precession formula 
(Eq. [17] with F = 0) and the exact calculation is obtained for 
Xt = 0.2 , but not for %r = 0.8 and 1.6, as expected from Eq. 36 
Figures [TO];b) and ^Uj^c) compare the exact calculation with 
the MF formula for A/^ = 50 and 100 spins. Excellent agree- 
ment between the MF formula and the exact calculation is ob- 
tained for xt = 0.2 and 0.8. Some differences are observed at 
Xt = 1.6. These differences decrease as N increases. 



1 2 3 4 5 



50 Particles 




Figure 10. vs initial polar angle 9i for an initial coherent spin 
state of N spins, t = IXarm in the legend is the total interaction period. 
The solid line is the exact calculation. The open circles are from the 
MF precession formula (Eq.[T7l with F = 0). (a) = 5. (b) = 50. 
(c)A^= 100. 
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